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Effects of turbulent pumping on stellar activity cycles 
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Stellar magnetic activity of solar like stars is thought to be due to an internal dynamo. While the Sun has been the subject 
of intense research for refining dynamo models, observations of magnetic cyclic activity in solar type stars have become 
more and more available, opening a new path to understand the underlying physics behind stellar cycles. For instance, 
it is key to understand how stellar rotation rate influences magnetic cycle period Pcyc. Recent numerical simulations of 
advection-dominated Babcock Leighton models have demonstrated that it is difficult to explain this observed trend given 
a) the strong influence of the cycle period to the meridional circulation amplitude and b) the fact that 3D models indicate 
that meridional flows become weaker as the rotation rate increases. In this paper, we introduce the turbulent pumping 
mechanism as another advective process capable also of transporting the magnetic fields. We found that this model is now 
able to reproduce the observations under the assumption that this effect increases as Q^. The turbulent pumping becomes 
indeed another major player able to circumvent the meridional circulation. However, for high rotation rates (SI ~ 5^q), 
its effects dominate those of the meridional circulation, entering a new class of regime dominated by the advection of 
turbulent pumping and thus leading to a cyclic activity qualitatively different from that of the sun. 
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, ^, 1 Introduction 



1.1 Magnetic cycles in Sun and Stars 

The Sun exhibits a magnetic period of about 22 years. Ob- 
servations based on the surface activity see n in Ca II- HK 
emission lines in the infrared and X-ray flux (' Saadl2009h re- 
vealed the existence of such cycles with comparable periods 
(from 5 to 25 years) in solar like stars, i.e. constituted of a 
deep convective envelope and a radiative interior. The avail- 
ability of more and more data on the magnetism of G&K 
stars provide us a new insight in understanding the global 
scale generation of their magnetic fields. This magnetism 
and regular activity is thought to be due to a magnetic dy- 
namo oper ating in the bulk and at the edges of the convec- 



tion zone jBaliunas et al. 1996' 
2004 ISaar & Brandenburg 1 999h 



Jouve et al. 2010a: Nandv 



In order to complement costly full 3D magnetohydrody- 
namical (MHD) simulations, a useful appr oach has been t o 
make use of the mean field dynamo theory iMoffatti (Il978h . 
This method has the advantage that it only deals with the 
large scale magnetic field, assuming some parametrization 
of the underlying small scale turbulence and magnetism. In 
these theories, the toroidal field is generally assumed to be 
generated at the base of the convection zone (BCZ) in a 
region called the tachocline, where both a radial and lat- 
itudinal shear act. What is more poorly known though is 
the source of the poloidal field to close the dynamo loop 
(i.e. Bpoi — ?> Btor — > Bpoi). Among the various existing 
mechanisms, one of the most promising one is the Bab- 
cock Leighton (BL) mechanism first proposed by iBabcock 



( Il96lb and further elaborated by iLeightonI (Il969h . In BL 
models, the poloidal field is generated at the surface by the 
transport and decay of tilted bipolar magnetic regions which 
are formed by twisted buoyant magnetic flux ropes devel- 
oped at the tachocline (IBrowning et al ]l2006l: iJouve & BrunI 
20091) . Synoptic magnetographic monitoring over solar cy- 
cles 21 and 22 has offered strong evidence for such dyn amo 



action toasi-Espuig et aI.l2010l:lWang et al.ll989[ll99lh . The 
BL flux transport dynamo models have thus been recently 
the favored one and have demonstrated some success at re- 
producing solar observations assuming either an advection 
(e.g. Dikpati & Charbonneau 1999) or a diffusion dominated 
(e.g. lYeates et al.ll2008l) regime for the transport of field 
from the surface down to the tachocline. Alternative a ef- 
fects such as helical turbulence may also be used in mean 
field dynamo models (ICharbonneatjl20 1 Oh but we will not 
do so here. 

In the past decades, dynamo action has been studied 
preferentially on the Sun, but some works have addressed 
the question of the applicability of such models to other so- 
lar type stars which possess different rotation rates a nd ac- 
tivity levels ( iBahunas et alJll996l: iNoves et al]|l984 . One 
difficulty of such observational programs is that they re- 
quire long term observations since stellar cycle periods are 
likely to be commensurate to the solar 1 1 -yr sunspots cy- 
cle period (or 22-yr for a full cycle including two polarity 
reversals of the global polo idal field). The biggest survey 
to date is from Mt. Wilson dSaar l2009h . but gathered data 
from the literatur e are now available f or solar type stars with 
enough statistics (iWright et al.ll201 ll) . The systematic anal- 
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ysis of these data revealed that the cycle period is shorter as 
the star rotates faster and that above a certain rotation rate 
the X-ra y luminosity satura tes (see below). The pioneering 



work of iNoves et al 

n 



(1984) found that Pc- 



P^with 



1.25 ± 0.5. However. [Brandenburg et alj (119981) ar- 
gued that an "Active" (with n = 0.8 ) and "Inactive" (with 
n = 1.15)branches, segregating respectively Young and old 
stars, could actually coexist. The Sun is found to actually 
lie in between. These scalings take into account the exis- 
tence of a primary (Hale) cycle and a secondary (Gleiss- 
berg) cycle. However, saturation of the X -ray luminosity 
limits this scaling to mod erate rotation rates dPizzolato et al.l 
iooi lWright et aLlEoT Jj). For G type stars this saturation is 
found for rotation rate above 35 km s^^, for K type stars 
at about 10 km and for M dwar fs around 3-4 km s^^ 
( Brownind2008 : Reiners et al .l2009h . ] low stellar magnetic 
flux scales with rotation rate is thus also important to un- 
derstand since it tells us how the magnetic field generated 
by dynamo action inside the stars emerges and imprints the 
stellar surface ( Rempel 2008) and if it actually saturates. 

In the framework of BL flux transport models including 
a meridional circulation (MC), the Pcyc — Prot relationship 
can be reproduced only if the meridional flow is propor- 
tional to the rotation rate of the star (Charbonneau & Saar 



2001 



Dikpatiet all 12001 



Nandyll2004t INandv & Martens 
2007[). However recent t heoretical work bv lBallot et al., C2007 ) 



and iBrown et alj ( 120081) indicate instead that the ampUtude 
of the meridional flow weakens as the rotation rate is in- 



creased. The recent work of iJouve et al.l (l2010al) (hereafter 
JBB2010) shows indeed that with such scaling, they cannot 
recover the observational trend of a shorter cycle period for 
faster rotation rates, unless a multicellular MC flow is as- 
sumed. We wish here instead to explore the influence of the 
turbulent pumping using a simple unicellular MC without 
invoking a more complex pattern. 



latter were able to reproduce rather well the large scale mag- 
netic field using only the alpha effect and differential rota- 
tion to drive the dynamo, the pumping effect was thought to 
be an unnecessary complication. 

A first approach showing t he importance of pumpin g 
in the solar cycle was made bv iBrandenburg etaL ( 1992 ). 
Since then, turbulent pumping has been a useful approach to 
tackle the problem of storage of magnetic fields. Indeed, for 
magnetic flux ropes to be buoyantly unstable and to emerge 
at the surface to form bipolar magnetic region with the ap- 
propriate tilt, numerical simulations h ave shown that their 
strength must be as high as Iff*— 10^ G dChoudhuri & Oilman 
1987; Fan 2004; Jouve & Brunll2007ah . One important lim- 
itation of this scenario is that lO'^ G represents a magnetic 
energy density higher by an order of magnitude than the ki- 
netic energy density, i.e. one needs to create super equipar- 
tition magnetic structures. Therefore, a stable layer is re- 
quired to store and amplify the magnetic fields. For this pro- 
cess to occur, differential rotation must be able to develop 
intense toroidal magnetic fields either within the tachocline 
or in the convection envelope. Pumping could be the trans- 
port needed to get the poloidal field down to the tachocline 
and could also maintain strongly buoyant structures from 
rising, thus helping them to become even stronger. 

More interestingly in the regard of our work, turbulent 
pumping h ave shown remarkable properties r egardi ng stel- 
lar cycles. Guerrero & de Oouveia Dal Find (l2008h (here- 
after GdG2008) have demonstrated that the magnetic cycle 
period Pcyc is no longer dominated by the meridional flow 
speed vq but instead by the radial turbulent pumping 7,. fol- 
lowing the relation : 



0.51^-0.05 



(1) 



1.2 Turbulent pumping in stellar dynamo 

Magnetic pumping refers to transport of magnetic fields in 
convective layers that does not result from bulk motion. One 
particular case is turbulent pumping. In inhomogeneous con- 
vection due to density stratification, convection cells take 
the form of broa d hot upflows surrounded by a network o f 
downflow lanes JCattaneo et al.lll99ll; iMiesch et al]|2008l) . 
In such radially asymmetric convection, numerical simula- 
tions show that the magnetic field is preferentially dragged 
downward (Tobi as et alJ200ll) . This effect has been demon- 



strated to operate in the bulk of the solar convection zone. 
A significant equatorward latitudinal component also arise 
when rotation becomes important, i.e. when the Rossby num- 
ber is less than unity. Turbulent pumping speeds of a few m 
can be reached acc ording to the numerical simulations of 



Kapvla et al 



(I2OO6'). Therefore, its effects are expected to 
be comparable to those of meridional circulation. 

In spite of those results, the effects of turbulent pumping 
rarely have been considered in mean field models. As the 



available in the range vo = [500; 3000], 7^0 = [20; 120], 
76»o = [60; 140] in cm s^*. Another interesting feature of 
this model is that the surface magnetic field no longer shows 
the strong concentration in the polar region that usually char- 
acterizes Babcock Leighton dynamo solutions operating in 
the advection dominated regime. This can be traced pri- 
marily to the efficient downward turbulent pumping that 
subducts the poloidal field as it is carried poleward by the 
meridional flow. 

With the STELEM code, we aim to study the effect of 
turbulent pumping on Pcyc and to see under which condi- 
tions the Pcyc — Plot relationship can be recovered. The 
paper is organised in the following manner In Section IH 
we describe the equations, the initial and boundary condi- 
tions, the ingredients of the model and the standard case 
in which the turbulent pumping is not included. Section [3] 
shows how the dynamo model behaves in the presence of 
turbulent pumping and we conclude in Section]!] 
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2 The model 

2.1 The model equation 

To model the stellar global dynamo operating in solar like 
stars, we start from the hydromagnetic induction equation, 
governing the evolution of the magnetic field B in response 
to advection by a flow field V acting against the magnetic 
dissipation characterised by the molecular magnetic diffu- 
sivity ijjn ■ 

— - Vx(VxB)-Vx {rj^V x B) (2) 

As we are working in the framework of mean field the- 
ory, we are interested in the large-scale magnetic field on 
time scales longer than the turbulent time scale. We ex- 
press both magnetic and velocity fields as a sum of a mean 
component (usually defined as a longitudinal average) and 
small-scale fluctuating component. For instance, the mag- 
netic field is decomposed as : 



B= (B) 



(3) 



Upon this separation and averaging procedure, the i 
duction equation for the mean component becomes 



m- 



9(B) 
dt 



= V X ((V) X (B)) + V X (v X b) 
- V X (77„V X (B)) 



(4) 



We drop now the averaging symbol () for the sake of 
clarity in the rest of the paper. A closure relation must then 
be used to express the mean electromotive force (emf) £ = 
(v X b) in terms of mean magnetic field, leading to a simpli- 
fied mean-field equation. If the mean magnetic field varies 
slowly in time and space, the emf can be represented in 
terms of B and its gradients 



Si — 



'ijk 



dxk 



(5) 



where and bijk are in the general case tensors con- 
taining the transport coefficients, and the dots indicate that 
higher order derivatives can be taken into account. Summa- 
tion over repeated indices is assumed. The tensors and 
bijk cannot, in general, be expressed from first principles 
due to the lack of a comprehensive theory of convective tur- 
bulence. In the kinematic regime where the magnetic energy 
is negligible in comparison to the kinetic energy, the most 
simple approximation is to neglect all correlations higher 
than second order in the fluctuations. This is the so-called 
first order smoot hing approxirnation (FOSA). In the pre- 
vious works (see ICharbonneaull20 1 01 for a recent review), 
they consider the simple case of isotropic turbulence, the 
tensor reduces into a single scalar giving rise to the a- 
effect. However, we consider in this paper the full tensor 
non-isotropic case for a. The emf can then in general be 
written as 



where a is a scalar referring to the standard a-effect. As 
we work with BL models, we assume instead the poloidal 
field to be generated at the surface, so that we will replace 
the aB term by an non local source term S (details are de- 
scribed below). The term 7 is the turbulent pumping and /3 is 
defined such that hijk — fieijk (with eijk is the Levi-Civita 
tensor). 

As 7 and /3 originates from the same velocity field, the y 
both in principle depend on v (see lKrause & Radleij[l980l) . 
In complex configurations as rotating convective spherical 
shells, it is not clear though how they are related in spite 
of numerous efforts in de termining precisely the turbulent 



transp ort coefficients (e.R. Brandenburg etal 



i2006t IKapyla et al.1 120091: lOssendrijver et al 



ml 



1990HCattaneo & Hughes 



200lh . In the 



models computed in this work, we thus do not consider any 
relationship between these two quantities. As can be seen 
in equations [TO] and [TT] (3 can be directly interpreted as 
an effective diffusion coefficient whose effects are known. 
Hence, we have focused on the effects of varying 7 for a 
fixed /3 even if a full parameter study in {7; /?} space should 
in principle be done. 

Working in spherical coordinates and under the assump- 
tion of axisymmetry, we write the total mean magnetic field 
B, velocity field V and the turbulent pumping (where we 
neglect 7^ because its amplitude is much less than the dif- 
ferential rotation) as : 



B(r, 6*, t) = V X {A^ir, 9, t)e^) + B^{r, 9, t)e^ (7) 
V(r,0) - Vp(r,0) +rsin0ri(r,^)e^, (8) 
7(r, 9) = 7p(r, 9) = 7r(r, ^)e, + 7e(r, e)ee (9) 

where Vp is the poloidal velocity field. Reintroducing this 
poloidal/toroidal decomposition of the field in the mean in- 
duction equation, we get two coupled partial differential 
equations, one involving the vector potential and the 
other concerning the toroidal field S^. The corresponding 
dimensionless equations are then : 



dA^ 
dt 



- f -J-]a^ + CsS{r, 9, B^) 
w 

+ — {C^rlrGr + C^eleGe) ■ V(n7A0) (10) 



dB^ 
dt 



,Brt 



1 d{voB^) d{r]/r,t) 



dr 



dr 



Z = (aB + 7 X B) - /3V X B 



(6) 



- i?etI7Vp • V(^) - i?eS^V • Vp 

\ w 

- B^V ■ {CyrlrGr + Cjejoee) 
+ Cntn(V X (zuA^e^)) ■ Vn 



(11) 
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where m ~ r sin 0,7] = rjm + j3, rjt is the turbulent 
magnetic diffusivity (diffusivity in the convective zone) and 
S{r,d,Bci,) the Babcock-Leighton surface source term for 
poloidal field (we neglect its contribution in the generation 
of the toroidal field compare to the shear applied by the dif- 
ferential rotation. We define the diffusive timescale as = 
Rq /rjt- Note that our velocity field is time-independent since 
we will not assume any fluctuations in time of the differen- 
tial rotation Q, or of the meridional circulation Vp . 

In order to write these equations in a dimensionless form, 
we choose as length scale the solar radius Rq and as timescale 
the diffusion time Rq/tji based on the envelope diffusivity 



with Tc = 0.7, di = 0.02, = 0.92 and C2 = 0.2. With 
this profile, the radial shear is maximal at the tachocline. 



Vt 



= 5 



This leads to the appearance of five 



control parameters Cq = Q,QR^/rit, Cs = sgRQ/rit and 
Re = VQRQ/?]t where nQ,so,vo are respectively the ro- 
tation rate at the equator and the maximal amplitude of the 
surface source term and of the meridional flow. For the solar 
rotation rate, we have f2o/27r = 456 nHz. We can also de- 
fine similar dimensionless numbers C^,. = JtoRq I Vt and 
C.yg = jeoRQ /vt for turbulent pumping. As we don't know 
how the Cy's are related, we allow for any relation between 
them. 

Equations [TOl andfTTI are solved with the STELEM code 
(Emonet & Charbonneau 1998, private communication) in 
an annular meridional plane with the colatitude 6 G [0, tt] 
and the radius (in dimensionless units) r e [0.6, 1] i.e. from 
slightly below the tachocline (r = 0.7) up to the solar sur- 
face. The STELEM code has been thoroughly tested and 
validated thanks to an internation al mean field dynamo bench 
mark involving 8 different codes dJouve et al I2OO8). At 6* = 
and 9 = TT boundaries, both A^f, and are set to 0. Both 

and are set to at r = 0.6. At the upper bound- 
ary, we smoothly match our solution to an external poten 



tial fiel d, i.e. we have vacuum for r > 1 (see Jouv e & Brun 



2007bl for more details). As initial conditions we are set- 
ting a confined dipolar field configuration, i.e. the poloidal 
field is set to sin 6'/r'^ in the convective zone and to below 
the tachocline whereas the toroidal field is set to every- 
where. All simulations have been carried with a resolution 
of 129x129 and the parameters used for the models are sum- 
marised in Table [T] 

2.2 The physical ingredients 

All quantities are in dimensionless values, and all the pro- 
files are normalised to unity such that the dimensionless pa- 
rameters set the strength of each ingredient. 

The rotation profile captures some realistic aspects of 
the Sun's angul ar velocity, ded uced from helioseismic in- 
versions (.Thompson et al.ll2003i) . assuming a solid rotation 
below 0.66 and a differential rotation above the interface (cf. 
Fig.ffll. 



1 



2 

X [1 - n 



1 + erf 



di 





Fig. 1 Left panel : Isorotation lines. Contours are reg- 
ularly spaced by 0.02 and spans from 0.8 to 1.0. Right 
panel : Meridional streamfunction. Contours are logarith- 
mically spaced by lO"'^. Positive (negative) values are in 
solid (dash) line. The thick dash line locates the tachocUne 
position. 



We assume that the diffusivity in the envelope rj is dom- 
inated by its turbulent contribution r/T in the convection 
zone whereas in the stable interior, rjc <C tjt- We smoothly 
match the two different constant values with an error func- 
tion which enables us to quickly and continuously transit 
from rjc to i]t in the vicinity of the core-envelope interface 
i.e. 



rjt 2 \ rjT 



1 + erf 



(13) 



with rjc = 5 lO^cm^s ^. 

In Babcock-Leighton (BL) flux-transport dynamo mod- 
els, the poloidal field owes its origin to the tilt of magnetic 
loops emerging at the solar surface. Thus, the source has to 
be confined to a thin layer just below the surface and since 
the process is fundamentally non-local, the source term de- 
pends on the variation of at the BCZ. A quenching term 
is introduced to prevent the magnetic energy from growing 
exponentially without bound (see Section |4|i. The expres- 
sion is then 



Sir,9,B^ 



1 + erf 
1 — erf 

1 



r - r2 
d2 
r- 1 



B^{rc,9,t) 

Bn 



(14) 



C2 C0S^{9)] 



(12) 



X cos a sm 



m9B^irc,9,t) 
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where r2 = 0.95 is the location of the max of the BL source 
term, d2 ~ 0.01 is the thickness of this layer and = 
lO^G. 



1 .000 r 



0.100 ^ 



'2008V In the series of models discussed in this paper, we re- 
strict ourselves to Babcock-Leighton flux-transport models 
that have a large single cell per hemisphere. As in the Sun, 
the meridional circulation are directed poleward at the sur- 
face and here they vanish at the bottom boundary (r = 0.6). 
This flow penetrates a little beneath the t achocline into the 



radiative interior as it is li kely to occur dBrun et al.ll201 1 



Garaud & Brummell 2008 ). To model the single cell merid 



ional circulation we consider a stream function with the fol- 



lowing expression from lDikpatil (1201 11) : 



0.010 ; 



r 



0.001 



\l-rb 



' — 9o) sin I TT 



1 - 



n 



(15) 



0.6 



0.7 



0.8 



0.9 



1 .0 - 



0.8 - 



0.6 - 



45 deg 
20 deg 
10 deg 



0.4 - 



0.2 



0.0 



where n = 0.65, = 0, T = 0.15, tq = 0.76, e = 
2.0, /3i = 0.316228, /Sa = 0.3. Vo is chosen such that the 
velocities are normalised to unity. We define p the density 
given by 

- -| - m 

p{r)= --0.97 (16) 
[r 

in which m = 1.5 The velocity components are derived 



through the relation pvp = V x 



0.6 



0.7 



0.8 
r/Ro 



0.9 



1 .0 



In our simula- 
tion, radial resolution spans from 0.6 to l.Oi?© correspond- 
ing to a density contrast of Ap ~ 120 

With this setup, the ratio between the maximal speed at 
the surface and the BCZ is Wsurf/i'BCZ — 109. This flow 
further vanishes at rf,. We chose this profile because it al- 
low us to control this ratio, which we found to be an impor- 
tant parameter for the ability of pumping to influence the 
magnetic period. The importance of pumping is enhanced 
as the MC amplitude is decreased at the BCZ, i.e. there is a 
stronger correlation between the magnetic cycle period and 
the amplitude of pumping. 

r,,, , J. 1. 1 J ^ £ -^u We turn now to the description of the turbulent pumping. 

The MC IS observed to be poleward at the sur face with a ^ i- i- o' 

J £ or, -1 . ji J JT^ o A ^- (he process we will focus on. In prmciple, its characteristics 

maximum speed of — 20 ms at midlatitude dBasu & Antial ^ . _j j- • „i • . _r ^, . 

l2010l) . but helioseismic inversion are currently unable to 



Fig. 2 Upper panel : Diffusivity profile as function of ra- 
dius. Lower panel : Babcock Leighton source term as func- 
tion of radius for different colatitudes. 



probe its amphtude in layers much below 0.95i?(7, as it pos- 
sesses we^CTfl2ws_^at are diffi cult to detect (Giles et al. 
1997l:lHaber et alJ2002bl:lHathaway & Rightmirel2010l) . Nev- 



can be determined by direct numerical integration of the 
equations of magnetohydrodynarnics (MHD). Se veral au- 
thors (e.g. [Brandenburg et alliggctfOssendriiver et alj200lh 



ertheless, this element is crucial in flux transport models in- 
corporating such process because it links the two sources 
of the magnetic field, namely the BCZ and the solar sur- 
face. In theoretical models, one usually assumes to conserve 
mass using a unicellular flow with an equatorward return 
flow at the BC Z. This has been re produced in recent so- 
lar simulations ( Miesch et al.ll2008l) when global flows are 
averaged over long intervals. However, other simulations 
tend to show that it is instead rather multicellular both in 
radial and latitud inal direction and highly time dependant 
jBrowning et al.l f2006: Brun & Toomre 2002). 

The nature of these flows in other stars is even less con- 
strained. In more rapidly rotating Suns, 3-D models presently 
indicate that the circulatio ns are likely to be multi-celled 
in both radius and latitude (IBallot et al.ll2007t iBrown et al 



have attacked this problem in order to calculate the dynamo 
coefficients, i.e. the a-effect and 7 -pumping. However, phys- 
ical conditions in the solar convection zone (e.g. high Reynolds 
number) prevent simulations from computing these coeffi- 
cients in an entire shell. For this reason, most of 3D MHD 
simulations are restricted to Cartesian boxes that represents 
only a small section (in both ra dius and lati tude) of stel- 
lar convection zone (see however ferown et"aL l 2010, for re- 
cent progress in computing transport coefficients in global 
model). The profile we use in this work a re inspired by one 
of the last estimation done bv lKapvla et al.. (2006.) . 



1 + erf 



1 - erf 



r~ 0.715 
0.015 
r-1.02 
0.05 
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79 



X COS^{0) 

25^5 
64 



(17) 



1 + erf 
erf 



0.74 



0.03 
r-0.94 



0.03 



X cos{0) sm^{9) 



(18) 



1 .0 



0.5 



0.0 



-0.5 



-1 .0 



30 45 45 
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90 



1 .0 



0.0 



-0.5 



-1 .0 - 



0.6 



0.7 



0.8 

r/Rc 



0.9 



Fig. 3 Profile of 79 (dash line) and 7^ (solid line) as func- 
tion of latitude at their maximal radial strength, i.e. at 1-0.85 
for both cases (upper panel), and as function of radius at 
their maximal latitudinal strength, i.e. at 9 = for 7^ and 
9 — acos(l/\/5) for7e (lower panel). 



2.3 The standard model 

We first compute a model without the pumping effect which 
we will refer as our standard case for direct comparison 
with forthcoming models. As we want to reproduce the Sun's 
properties, we have chosen our parameters for case S (see 
Table [1]) such as to get the correct cycle period. We choose 
Cs such that it is well above the threshold for dynamo ac- 
tion (C"'' = 9.6), enabling us to see the fields in the regime 
of a well established flux transport dynamo. 

We show in Fig.|4]the radial field at the surface and the 
toroidal field at the BCZ, i.e. at r = O.7i?0. It can be di- 
rectly compared to the solar butterfly diagram if we iden- 
tify the toroidal field at the BCZ as the source of buoyant 



magnetic flux tubes rising radially through the convection 
up to the surface seen as active regions. For this, the ris- 
ing time must be very short compare to the magnetic cycle. 
But according to models, the rising time spans from months 
to dozens of years (see review of Fan.2009,). In particular. 



Jouve et al.l (120 1 Obi) found that even a short delay is impor- 



tant for a modulation in the cycle but not on the cycle period 
itself. We will thus not address this problem in this work. 



0) 




100 



120 



180 



200 



140 160 
Time [years] 

Fig. 4 Butterfly diagram, i.e. time-latitude slice of the 
toroidal field at the BCZ and of the radial field at the sur- 
face. Red (blue) colors indicate positive (negative) magnetic 
fields. Contours are logarithmically spaced with 2 contours 
covering a decade in field strength. The color range spans 
from Max to Min of the magnetic field strength. The vertical 
dashed line corresponds to the epoch of reversal of toroidal 
field, the plain line correspond to the epoch of reversal of 
poloidal field at the poles from negative to positive polarity 
and the dash-dotted line coiTesponds to the positive maxi- 
mum of toroidal field near the equator. 



With this typical model, we are able to reproduce sev- 
eral aspects of the solar cycle, notably its period of approx- 
imately 22 years, a strong equatorward branch for toroidal 
field restricted to low latitudes, a phase shift of 7r/2 between 
the surface polar field and the deep toroidal field, so that 
the polar field changes its polarity from negative to positive 
when the toroidal field is positive and maximal in intensity 
near the equator. 

As described in Section l2!2l the amplitude of MC is very 
low at the BCZ compare to the velocity at the surface. In the 
advection dominated regime, we found that this strongly af- 
fects the magnetic cycle periods because it is the dominant 
mechanism capable of transporti ng the toroidal field toward 
the equator. The recent work from lPipin & Kosovic hev ( 2 01 lit 
confirm the great importance of the speed and the depth of 
the return flow. To compensate for this effect and to recover 
the 22-year magnetic cycle, we choose a maximal velocity 
at the surface (^ 32 m s^^) slightly higher than what is usu- 
ally assumed in simulations (^ 20 m s^^). However, merid- 
ional flow is a very time-dependant process w hich can reach 
values as high as what is used in this work (IBasu & Antia 
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2010h . Also, the strong equatorward branch for the toroidal 
field is the signature of the drag of the toroidal field by equa- 
torward MC at the BCZ and thus clearly shows the domi- 
nating effect of field advection over diffusion which can be 
seen in the scaling law : 



P< 



eye 



0.05 „0.07„, -0.83 



(19) 



in the range So = [4.7; 100] cm s~\ = [1200; 4500] cm 
and Q.0 = [1.4 x 10-*^; 7.2 x 10"^] Hz. We directly see that 
the meridional circulation speed dominates over Cs and Cq. 
As dominates the scaling, it is crucial to assess how 
scales with f^n- We rely on 3D numerical simulation from 



Table 1 Summary of the control parameters defined in 
Section 12.21 for the different cases studied in this work. S 
and R stands for Standard and Reference, i.e. they include 
(respectively dot not include) turbulent pumping. The trail- 
ing number in the case name gives the rotation rate Q, in 
units of r^Q. The S and R cases are at the solar rotation rate. 
The advective control parameters (R^, C^r and C^e are nor- 
1 malized such that their relative strength can be directly esti- 
mated. 



Brown et alj (120081) which studies the influence of rotation 
rates on solar like stars. They found that MC decreases with 
the rotation rate as vq (x Qq^'^^. This is not intuitive as one 
could expect that the meridional circulation increases with 
the rotation rate. A careful study of the vorticity equation 
shows that it actually weakens with rotation rate as more and 
more kinetic energy is being transferred to longitudinal mo- 
tions at the expense of meridional kinetic energy. Assum- 
ing this relation in our simulations, we obtain naturally that 
Pcye « f^o-^i in the range f^o = [l.Ox 10^^; 1.5x 10^^] Hz. 
This is in agreement with what has been found in JBB2010. 

As stated above, the current model reproduce quite well 
a solar-like butterfly diagram, but is still unable to repro- 
duce the Plot — Peyc relationship because of the strong in- 
fluence of MC in these models, especially at the BCZ. To 
address this issue, JBB2010 tried to incorporate a multicel- 
lular MC as observed i n 3D numerical sirn ulations of e.g. 



BaUot et al.l (l2007h and iBrown et all (|2008|) . and were in- 
deed able to recover the observational trend. In this work, 
we propose another solution, that is to keep a unicellular 
MC but to introduce the turbulent pumping as a new mech- 
anism to transport the poloidal field from the surface down 
to the tachocline. GdG2008 offered insights on the possibil- 
ities of such process for a shallow MC (see Section [L2] i. We 
wish here to verify if their results holds for a deeper MC and 
under which conditions the turbulent pumping can shorten 
the advection path driven by the meridional circulation. 

3 Influence of turbulent pumping on stellar 
cycles 

3.1 Reference case 

We turn now to the models with turbulent pumping. As be- 
fore for the standard case, we define here the reference case 
R corresponding to the Sun, i.e. at the solar rotation rate 
(cf Fig. |5]l and reproducing the solar features described in 



Section 12.31 We rely on the estimation from iKapyla et al 

for the pumping amplitude with jro = 40 cm s 
and 7eo = 100 cm s^^. We sum up the physical parame- 
ters for the reference case (i.e. with turbulent pumping) at 1 
solar rotation rate as R in Table [T] 

With the profiles described in Section 12.21 the latitudi- 
nal pumping component increases the total advective speed 



Case 






IW 


IW 


-Pcyc 

(years) 


-Dpol/i'tor 


S 


2.78 


4.50 








22.5 


2.57 


R0.7 


2.00 


2.67 


0.0289 


0.0723 


26.4 


1.51 


R0.9 


2.50 


2.41 


0.0452 


0.113 


24.2 


0.999 


R 


2.78 


2.30 


0.0557 


0.139 


22.7 


0.701 


Rl.l 


3.00 


2.22 


0.0650 


0.163 


21.2 


0.536 


R1.3 


3.50 


2.07 


0.0885 


0.221 


18.7 


0.328 


R1.4 


4.00 


1.95 


0.116 


0.289 


16.3 


0.228 


R1.6 


4.50 


1.85 


0.146 


0.366 


14.5 


0.174 


R1.8 


5.00 


1.76 


0.181 


0.452 


13.0 


0.148 


R2.0 


5.50 


1.69 


0.219 


0.547 


11.7 


0.127 


R2.5 


7.00 


1.52 


0.354 


0.885 


9.22 


0.104 


R2.9 


8.00 


1.43 


0.462 


1.16 


8.22 


0.100 


R3.6 


10.0 


1.29 


0.723 


1.81 


6.76 


0.0706 


R4.3 


12.0 


1.19 


1.04 


2.60 


5.14 


0.0539 


R5.0 


14.0 


1.11 


1.42 


3.54 


4.51 


0.0377 



resulting in a lower MC amplitude (^ 17 m s~^) to keep 
a 22 year period, in a better agreement with temporally av- 
eraged observations. The strong equatorward branch of the 
toroidal field appears at slightly higher latitudes 65° with 
the regions of strongest magnetic intensity confined in a 
smaller area near the poles. We have the opposite situation 
at the surface where the effective speed is reduced expand- 
ing the poloidal field to lower latitudes. At the poles, the ra- 
dial pumping drags the poloidal field down to the tachocline 
and thus leads to a lower concentration of surface magnetic 
fields at the poles. This effect is in our deep MC not as strik- 
ing as what has been reported by GdG2008. The phase re- 
lation of 7r/2 between the polar and the toroidal fields de- 
scribed in Section 123] are well reproduced in this case (see 
Fig.E). 

Figure |6] shows the behaviour of the magnetic field in 
this model. Here the magnetic field follows the advective 
path created by the meridional flow and consequently this 
large scale flow plays a key role in these dynamo solutions. 
In this reference case at I^Iq, pumping amplitudes are not 
sufficiently high to advect the magnetic fields away from 
MC flow. 



3.2 Varying the rotation rates 

We wonder now under which conditions the turbulent pump- 
ing can solve our initial problem, namely the scaling rela- 
tionship between the magnetic cycle period Pcye with the 



rotation period P^ot- A least square fit through our data gives 

(20) 



Pcye OC V, 
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Fig. 5 Butterfly diagram for the reference case. Same 
color code as in Fig.|4] 




Time: 0.65 




Fig. 6 Evolution of the vector potential (black lines) and 
of the toroidal field (color contours) during half a mag- 
netic cycle. Red (blue) colours indicate positive (negative) 
toroidal field and plain (dotted) lines indicate clockwise (an- 
ticlockwise) poloidal field lines. Time is expressed in years. 



in the range of wq = [1000; 4500] cm s-\ 7^0 = [20; 100] 
cm s^^, 7eo = [50; 300] cm s~^. This result is not in com- 
plete agreement with the work of GdG2008. We found that 
indeed the turbulent pumping becomes a major player in 
setting the magnetic period, but the quantification of its in- 
fluence remains different. First, the MC is still the dominant 
effect and the radial pumping component is not as important 
as in GdG2008. Second, the effect of 761 is not negligible. 
This supports the idea that the latitudinal advection process, 
and especially at the BCZ, is an important ingredient in ad- 
vection dominated BL models, capable of transporting the 
toroidal magnetic field from the pole toward the equator. 
This difference may come from their choice of a shallow 
MC with almost zero velocity at the BCZ. 

A simple look at the scaling law|20]gives that if we want 
to recover the observational trend (again, we assume that 
vq (X ri~°-^^), and assuming that 7r/7e remains constant, 
the pumping effect should roughly scales as f2§. A first es- 



timate done in lTobias et al.l (120011) have shown that the tur- 
bul ent pumping actuall y decreases with rotation rate. Later 



on. 



Kapvla et al. I (l2009h found on the contrary that the rota- 



tion rate have almost no effect on the vertical pumping. This 
lets the scaling as an open question. In order to verify if this 
scaling holds in solar type stars, work are currently done in 
3D MHD simulation jMatt et al.ll20Tll) . 



30 p 
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3 
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Fig. 7 Magnetic cycle period as function of the rotation 
rate in models including turbulent pumping. Solid line is a 
least square fit of the simulated data. 



Under these assumptions, we found indeed that Pcyc oc 
in agreement with the observations as expected. On 
Fig. Ill we present the rotation rates ranging from 0.751© up 
to Sri© with a least square fit of the data. Our result does not 
hold outside this range because a systematic period does not 
emerge for 51 > 551©. For fl < O.Til©, the pumping is neg- 
ligible and lets MC imposing that Pcyc increases with ft. We 
are thus back to a strong dependency of Pcyc with MC am- 
plitude. We show in Fig. [8] 3 representative cases at 0.7il©, 
2flQ and SJl© . We see that the equatorward branch becomes 
shorter and shorter as the rotation rate, and so the pumping, 
is increased. Also, the surface magnetic field becomes ho- 
mogeneously distributed in latitude thanks to the increase 
of both components of the turbulent pumping. However, the 
determination of a period for the most rapidly rotating stars 
(fl > 4110) becomes difficult as the butterfly diagram is 
affected by more and more intermittency and small scale 
structures as can be seen in the bottom panel of Fig.[8] This 
explains the last points of Fig. |7] to lie sUghtly away from 
the trend. 

One can also note the appearance of a modulation of 
both cycle strength and period as soon as 2110 as seen in the 
middle panel of Fig. [8] We might not compare this with long 
term modulation known as Gleissberg cycles. For 2D mod- 
els in the kinematic regime (as in this work), the Lorentz 
force exerted by the magnetic fiel ds on the velocity field, th e 
so-called Malkus-Proctor effect (iMalkus & ProctojUQTSb . 
is not included. Such a feedback from the large scale Lorentz 
force has been shown to have strong effect on for instance 
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Fig. 8 Butterfly diagram for 3 representative cases : 
O.TrJ©, 2.00© and 5ilQ. Same colorcode as in Fig.H] Afl of 
the figures share the same color scale : between —5 x lO'^G 
and 5 x lO^G for and between -9 x lO^G and 9 x lO^G 
for B^. We show the phase relations for the top and middle 
panel only as it is not possible to define them for the highest 
rotation case. 



increasing the rotation rate. The solar butterfly diagram fea- 
tures almost vanish in very small structures. The poleward 
and equatorward branches barely exist and the magnetic 
field is completely homogeneously distributed over all lati- 
tudes and do not present a strong concentration of poloidal 
fields at the poles. In these cases, the advection is no longer 
dominated by MC but by pumping, suggering that we are 
entering in a new class of regime for BL flux transport mod- 
els. 

We turn now to the ratio between the maximal value of 
the poloidal field at the surface and the maximal value of 
the toroidal field at the BCZ, Bpoi/Stor- This ratio is found 
to decrease with the rotation rate (-Bpoi/-Stor oc as 
seen on Fig.|9] This is in reasonable agreement with the ob 



servations of rapidly rotating solar like stars by iPetit et al, 



( 120081) where rapid rotators host a large scale toroidal com- 
ponent in their surface field whereas the magnetic field is 
mostly poloidal for low rotation rates. In the range of avail- 
able observations, we found systematically a lower ratio. 
This is not surprising as we estimated iJtor at the BCZ where 
it is generated, and hence where it is the strongest whereas 
observers have access only at the weaker surface toroidal 
field. To evaluate the surface toroidal magnetic field in our 
simulations, a first approach would be to take the value close 
to the surface. A typical ratio between the BCZ (r — 0.7) 
and the surface (r = 0.98) toroidal fieldfor the reference 
case is found to be Bf^'^/Bl'^^ ~ 60. 
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Fig. 9 Bpo\/Btoi ratio as function of the rotation rate. 
Solid line is a least square fit of the data. 



torsional oscillations kempefcOOel) . l ong term modulation 
JPhillips et al.l2002h and intermittency dMoss & Brookel2000h 
In the current simulations instead, the modulations are due 
to the pumping mechanism which advects the magnetic fields 
against the MC flow. Thus, the phase becomes more and 
more modified between Bpoi and i?tor up to a point where 
ephemeral dynamo loops appears, generating multiple pe- 
riodicities. It becomes even more dramatic when we keep 



4 Conclusion and perspectives 

The aim of this study was to address the current issue en- 
countered by 2D mean field models in which the magnetic 
period dec reases with the rotation rate, co ntrary to the ob- 
servations (ISaaij|2009t IWright et al.ll201 11) . Earlier work of 
JBB2010 showed that this behaviour was due to the large in- 
fluence of the MC. More spec ifically, they have us ed the re- 
sults of recent 3D simulations iBrown et al.l ( 2008 ) in which 
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meridional circulation amplitude decrease with rotation rate, 
and came to the conclusion that the observed Pcyc — f rot re- 
lationship cannot be reproduced, unless a multicellular MC 
is considered. 

The idea of this work was to tackle this problem un- 
der another angle introducing the turbulent pumping mech- 
anism which has been shown to have interesting proper- 
ties on the magnetic period (GdG2008). We have then per- 
formed 2D BL flux transport simulations with the STELEM 
code in the advection dominated regime. We first computed 
a standard model without turbulent pumping but producing 
solar characteristics. We found that in such models the mag- 
netic period is indeed very sensitive to the MC amplitude 



0.05 ^0.07„, -0.83 



), and is therefore not able to re- 



cyc '^'0 ■'0 '^O 

produce the observations, confirming the results of JBB2010. 

In the presence of turbulent pumping however, the MC 
is no longer the only process capable of influencing -Pcyc- 
We found that the latitudinal speed at the tachocline drives 
the magnetic cycle period. The weaker the MC, the stronger 
the influence of pumping on the magnetic cycle period. A 
reasonable equatorward flow at the BCZ 0.2 ni s^^) with 
a velocity contrast between the surface and the tachocline of 
- 109 gives the scaling Pcyc oc V" '*"7rfi°'^°7eo°'^^- The 
observational trend can be thus recovered only if 7 oc 
which has not been reported yet in previous C artesian simu- 
lations dKapvla et al.l2009HTobias et al.l200lb . However 3D 
MHD simulations in full spherical sh ells are currently under 
development to test this assumption ( Matt et aklbOl 1 ). 

Another successful feature resulting from the presence 
of turbulent pumping is their ability to reduce the strong 
concentration at the poles. On one hand, the strong radial 
component drags the surface field down to the tachocline. 
On the other hand, the latitudinal component is equator- 
ward everywhere (and therefore at the surface) expanding 
the strong poloidal field to lower latitudes. As we increase 
the rotation rates, the pumping becomes stronger and this 
dilution is enhanced until the surface field becomes com- 
pletely homogeneous. However, at this level of rotation, the 
butterfly diagram is strongly affected by the pumping which 
has become the dominant advective process with a greater 
amplitude than the MC. For instance, it does not show strong 
equatorward nor poleward branches. We enter here in a dif- 
ferent dynamo regime, in which we are dominated , not by 
MC, but by turbulent pumping, with characteristics different 
than the Sun. 

Although this new ingredient looks promising with re- 
spect to our study, turbulent pumping amplitude must be as 
high as few ^ 10 m s^^ for ^ Sil©. Such high values have 
not been reported yet in direct numerical simulations. A way 
out would be to mix the two proposed approaches by con- 
sidering in the same time both pumping and multicellular 
MC. The latter naturally arised in 30 MHD global simula- 



tions jBaUot et al.l2007tlBrown et alkOOSllBrun & Toomre 



Finally, long term modulations have not been consid- 
ered in this work since no nonlinear effects beside quench- 
ing has been used. We intend to take into account those ef- 
fects such as Malkus Proctor term in a future work since 



2002h in which they last for several years and are supported 
by obs ervational evidences via local helioseismology tech- 
niques dHaber et alJl2002ah . 



they are potentially important ( iMoss & Brookd 12000,) . We 
also intend to extend the resolution domain by including an 
atmosphere which would give a prediction for the external 
toroidal field, but also a bottom bo undary condition for stel- 
lar wind models (iPinto et al 

Studying stellar magnetism with mean field models is 
instructive as we can extract information on the sensitivity 
of magnetic cycles to parameter change, such a study being 
very comput ationally expensive and delicate with 3D MHD 
simulations (IBrown et al.l 1201 lb . However, in the light of 
this work, a subtle combination of different processes must 
be considered to account for the many aspects of the obser- 
vations. More accurate data would be of great help to un- 
derstand the complex underlying physics of stellar dynamo. 
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